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We have searched for the supersymmctric partner of tlie top quark (stop) in 107 pb^^of pp colh- 
sions at y/s= 1.8 TeV collected by the Colhdor Detector at Fermilab (CDF). Within the framework 
of the Minimal Supersymmetric extension of the Standard Model (MSSM) each of the pair-produced 
stops is assumed to decay into a lepton, bottom quark and supersymmetric neutrino. Such a sce- 
nario would give rise to events with two loptons, two hadronic jets, and a substantial imbalance of 
transverse energy. No evidence of such a stop signal has been found. We calculate a 95% confidence 
level (C.L.) upper limit on the stop production cross section, which excludes stop masses in the 
region (80 < < 135 GeV /(?) in the mass plane of stop versus sneutrino. 
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One of the most promising theories beyond the Stan- 
dard Model (SM) [1] is the Minimal Supersymmetric 
Standard Model (MSSM) [2]. It predicts that each SM 
particle has a superpartner (sparticle) with the same 
quantum numbers, except for spin which differs by one 
half unit. Experimental results indicate that supersym- 
metric (SUSY) particles are generally not as light as their 
SM partners. SUSY, therefore, is broken at or above the 
electroweak scale, and we treat the sparticle masses as 
free parameters. Due to the large top quark mass, there 
may be a large mixing between the superpartners of the 
left and right helicity states of the top quark [3]. This 
can lead to substantial mass splitting of the stop mass 
eigenstates {ii, ^2) with the lighter one (denoted i from 
now on) potentially being the lightest squark. 

Stop-antistop pairs {tt) arc strongly produced in the 
pp collisions at the Fermilab Tevatron if kinematically 
accessible. The production cross section has been cal- 
culated using QCD in the next-to-leading order (NLO) 
approximation [4]. For a given stop mass (wj) the cross 
section depends only weakly on the other parameters of 
the MSSM. In the mass region of interest to our search 
(m^ =80- 140 GeV/c^), the cross section drops from 44 pb 
to 1 pb. 

We assume SUSY i?-parity [5] conservation, from 
which the stability of the lightest supersymmetric parti- 
cle (LSP) follows. All SUSY particles, including the stop, 
eventually decay into this LSP. Stop decays into the top 
quark are kinematically not accessible in our region of 
interest due to the high top mass (m^ < mt). For the 
stop decay into a bottom quark and an on-shcU chargino 
(xf ), only a very small window of opportunity remains at 
the Tevatron due to the high xf mass limit from LEP2 
[6]. Another possible 2-body stop decay would be the 
flavor-changing, t cxi, decay [7]. It would proceed via 
higher order loop diagrams and is thus highly suppressed. 
The 3-body decay into a charged supersymmetric lepton, 
t lub, is closed for most of the stop region currently 
within the reach of CDF because of the slepton mass 
limit of LEP2 [6]. The existing mass limit of the su- 
persymmetric neutrino, ma > 45 GeV/c^ [8], leaves the 
decay into sneutrino, i — > IVb, open. This 3-body decay 
proceeds via a virtual chargino, and is expected to yield 
equal e, /x, and r branching ratios. 

Stop pair production with the i — > lub decay will yield 
two leptons with opposite electric charge, two hadronic 
jets from the bottom quarks and considerable transverse 
energy imbalance (EIt) in the detector [9] due to the es- 
caping sneutrinos. CDF has reported earlier on an anal- 
ysis based on B identification [10]. In this Letter we use 
dilepton events. Only a few SM processes yield dileptons 
and can thus potentially mimic our stop signature. The 
most significant ones are (1) tt production with leptons 
from W and/or bottom decays; (2) heavy flavor, i.e. bb 
and cc with semileptonic decays; (3) Drell-Yan produc- 



tion with hadronic jets from higher order processes; (4) 
diboson production, WW, WZ and ZZ; (5) lepton pairs 
from the decay of vector mesons, such as J/V' and T; 
(6) events without two genuine prompt leptons, where a 
hadron is misidentified as a lepton, or decays in flight to 
a lepton. 

The search presented here is based on 107 pb~^of pp 
collisions at ^/s = 1.8 TeV collected by the Collider De- 
tector at Fermilab (CDF) during the 1992 to 1995 run- 
ning period of the Tevatron. A detailed description of the 
CDF detector can be found in Ref. [11]. Online triggers 
selected approximately 6.4 million single lepton events 
and an additional 3.3 million dilepton events. All of those 
events have been reconstructed, and 13,295 events were 
selected as a dilepton sample, by requiring at least one 
tight electron (E^ > 10 GeV, jr/j < 1.0) or muon (pT > 
10 GeV/c, jr/l < 0.6) candidate, and a second loose elec- 
tron (Et > 6 GeV, [ryl < 1.0) or muon (pt > 6 GeV/c, 
\ri\ < 1.0) candidate. No explicit tau lepton identifica- 
tion was done, but taus can enter the search sample if 
they decay leptonically. Electrons are identified by en- 
ergy deposition in the electromagnetic calorimeter with a 
track of corresponding energy in the central drift cham- 
ber (CTC) pointing to it. Muons are identified by track 
segments in both the CTC and the muon drift chambers 
that are located behind 4.5 to 10 interaction lengths of 
absorber. Standard lepton identification cuts are used 
and described elsewhere [12]. Each lepton is required to 
be isolated, i.e. we require the total px of all other tracks 
within a cone AR = ^(A7?)2 + (A(/))2 < 0.4 around the 
lepton's track not to exceed 4 GcV/c. The jets were 
reconstructed with a cone algorithm with cone radius 
AR = 0.7 [13]. We require at least one jet in the central 
region of the calorimeter (|?7| < 1.0) with Et > 15 GeV, 
that is separated by AR > 0.7 from both leptons in the 
event. For increased efficiency, we require only one of the 
two jets to be identified. Sequential B decays. J/ih, T 
and Z events were removed reqiiiring the invariant dilep- 
ton masses mw > 6 GeV/c^ or mu > 12 GeV/c^ and ex- 
cluding 76 < niii < 106 GeV/c^ (where prime indicates 
any mixture of e and /i flavors and no prime indicates 
same-flavor dileptons). At the preselection level we start 
with > 15 GeV. Experimental backgrounds, like elec- 
trons from conversions and muons from cosmic rays are 
removed with additional cuts [14]. 176 events fulfill the 
above preselection requirements. 

To estimate the number of SM and stop events in the 
sample, events of the various physics processes are gen- 
erated by ISAJET [15] and simulated for the CDF detec- 
tor. We have used CTEQ-3 parton distribution functions 
(PDF) [16]. The stop production cross section was cal- 
culated with PROSPINO [17] and the ISAJET cross sec- 
tion was adjusted accordingly. We have generated events 
over a large range of stop (80-140 GeV/c^) and sneutrino 
(45-90 GeV/c2) masses. 

The Drell-Yan and ti production cross sections were 
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normalized to CDF measurements [18]. The Monte Carlo 
(MC) bb and cc cross sections were verified by inclusive 
electron-muon samples [19]. The B^B^ oscillation effect 
was added based on the CDF measured inclusive x mix- 
ing fraction [20]. The diboson production cross sections 
of the MC were scaled to those of NLO calculations [21]. 

For low pt leptons the contribution due to misidcn- 
tifi-cation can be significant and is calculated in two 
steps [14]. First we measure in various data samples 
the so-called "fake lepton probabilities" (momentum- 
dependent, separately for electrons and muons, and de- 
pendent on detector region). These "fake lepton proba- 
bilities" include hadrons being misidentified as electrons 
or muons, and also include leptons from in-flight decays 
of pions and kaons. The data samples used were (i) a 
minimum-bias trigger sample and (ii) a sample with a 
50 GeV jet trigger threshold. We remove events in those 
samples similarly to the dilepton analysis, and, to be un- 
biased, we exclude all hadrons associated with a jet re- 
quired by the online trigger. We measure misidentifica- 
tion probabilities between 0.4% and 7% for both e's and 
/i's [14], and the results are consistent between the two 
data samples. 

Second, in a single lepton sample we use these "fake 
lepton probabilities" successively on each other track 
in the event to simulate dilepton events. We use the 
"fake lepton probabilities" to simulate both the number 
of misidentified-lepton events as well as their kinematic 
properties. 

The major background to the preselection sample 
comes from heavy flavor production, with about a quarter 
of the events having leptons of the same charge. Another 
significant background comes from Drell-Yan processes. 
In those events the comes from r decays or jet and 
lepton energy mismeasurements due to uninstrumented 
detector regions. We expect a total background of 155 
± 55 events, while a stop and sneutrino mass combi- 
nation of 100 and 75 GeV/c^ would contribute 24 ± 9 
events. Table I shows the expected contributions in de- 
tail for like sign and opposite sign leptons. To verify our 
background calculation further, we compare kinematic 
distributions of the data and the expected background. 
Figure 1 shows a few such distributions. Top and di- 
boson production yield generally more energetic leptons 
than 66, cc or misidentified leptons. The px distribu- 
tions of the leptons show that both high and low px 
lepton sources agree well with the observed data. The 
observed FSt distribution agrees both at low Fix, where 
detector effects dominate, and at high Fix, where neutri- 
nos from W and Z bosons determine the spectrum. We 
also note that the parton shower MC describes well the 
observed jet multiplicity. From the signal to background 
ratio, it is clear that the preselection sample does not 
have sufficient sensitivity to answer the question of stop 
pair production. In contrast to an earlier search [22] we 
select a kinematic region in which we expect higher S/B 



where S is the stop signal and B is the background which 
includes SM processes and misidentified leptons. This 
decreases the systematic uncertainty associated with low 
Pt leptons in our analysis. 

TABLE I. Data, expected backgrounds for the prese- 
lection sample, and expected stop signal for mi {ma) = 100 
(75) GeV/c^. The stop event acceptance is 2.5% at this stage. 



Source 


OS 


LS 


Drell-Yan 


52.2 ± 13.7 


0.4 ± 0.4 


bb, cc 


43.5 ± 32.1 


16.4 ± 17.6 


tt 


9.5 ± 2.9 


0.6 ± 0.2 


WW, WZ, ZZ 


3.8 ± 0.9 


0.4 ± 0.1 


Misidentified Leptons 


16.3 ± 4.4 


12.4 ± 3.4 


Total Background 


125.2 ± 46.7 


30.1 ± 18.4 


Data 


128 


48 


Expected tt 


22.6 ± 8.9 


1.0 ± 0.4 
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FIG. 1. Data and expected background after preselection. 
Tight and second lepton transverse energies, missing trans- 
verse energy (for comparison we also show the missing Et 
distribution for a 10 times stop signal of mj (mj) = 100 
(75) GeV/c^) and jet multiplicity distributions shown for 
events with opposite charge leptons. The last high bins con- 
tain overflows. 



In less than 5% of stop events the two leptons will have 
the same electric charge due to the semileptonic decay of 
one of the 5-quarks. However, 20% of the SM background 
yields like-sign (LS) lepton events. We thus focus our 
search on events with opposite-sign charge (OS) leptons. 
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For i?p-conserving supcrsymmctry wc expect large miss- 
ing energy from the rather heavy sneutrinos. In Fig. 1 
we see most of the background events clustering at low 
missing E^. A cut of 30 GeV removes 77% of the 
SM background but keeps about 57% of the stop events. 
Energy mismeasurement of leptons, or the presence of 
neutrinos from Drcll-Yan t decays, would cause the lep- 
tons (and the dilepton system as well) to be aligned with 
the E!t direction. This is not typical for the signal, where 
we expect true from the sneutrinos in the stop de- 
cay. We therefore require the azimuthal angles between 
E!t and the individual leptons and the dilepton system 
A(Af^, A(^f^, and AtA^T to be larger than 30°. 

In Drell-Yan plus jets events or when bb or cc events 
originate from gluon splitting (initial or final state) 
events, the two leptons balance the jets in the transverse 
plane. We veto events where the angle between either lep- 
ton and the most energetic central jet, and A0j^\ 

is larger than 90°. 

Events from top pair production pass the above cuts 
with efficiencies similar to stop pair events and are now 
the dominant source of SM background. In top events 
the leptons come from W decay and are very energetic. 
In the case of stop, we have 3-body decays containing a 
very heavy sneutrino. The amount of available energy in 
the decay depends on the stop-sneutrino mass difference, 
Amj_,j. For small mass difference, the leptons and jets 
arc quite soft and a large fraction of the event energy es- 
capes detection through the sneutrinos, unlike a ti event. 
For best stop sensitivity at small Amf_p we require the 
scalar sum of lepton Pt, Pt + Pt — GeV/c, and the 
Pt of the dilepton system, prji ^ < 30 GeV/c. Although a 
large amount of energy escapes undetected, the sneutri- 
nos tend to be back-to-back, thus reducing the measured 
t-T We also require the sum of the most energetic central 
jet Et and the missing Et, EJ""^ + E^t < 160 GeV. 

For large stop-sneutrino mass difference, the leptons 
are more energetic and we can increase our lepton px 
requirement to 10 GeV/c without much loss in stop ef- 
ficiency. However, leptons and jets are still significantly 
softer than in tt events. Wc place the same jet, miss- 
ing Et, and lepton requirements as at small Amf_p, 
EJ^* + ESt < 160 GeV and p^^ + pf^ < 75 GeV/c, but 
loosen the requirement on the dilepton px to p?^'^ < 
55 GeV/c. 

Table II shows the expected number of stop events 
for the two search regions. We start our search at stop 
masses of 80 GeV/c^ to overlap with previous LEP limits. 
Near the kinematic limit of the stop decay, = mj>+m{,, 
lepton and jet energies become very soft, limiting our 
stop detection capabilities. At high stop mass our sen- 
sitivity is limited by the steeply falling production cross 
section. In the region of interest to this search the final 
stop event acceptance varies between 0.3% and 2.3%. 

The biggest source of uncertainty on the number of 



expected stop events arises from the choice of the rcnor- 
malization and factorization scale, Q^, which character- 
izes the amount of energy transferred during the colli- 
sion. The is reduced (due to the sneutrinos being 
more back-to-back) when Q is increased, and the jet Et 
gets softer when Q is decreased. By varying Q by a fac- 
tor of two up and down, we determine the uncertainty 
due to the choice of to be 32%. Other significant 
sources of uncertainty are: (1) the choice of PDF (11%) 
(2) the absolute energy scale of the detector (11%) (3) 
the amount of gluon radiation (7%) (4) trigger, lepton 
and isolation efficiency (5%), and (5) the luminosity mea- 
surement (4%). The statistical uncertainties of the MC 
samples are about 8%. Combining the statistical and sys- 
tematic uncertainties we obtained a total uncertainty of 
38% for the signal expectation. Similarly, wc evaluated 
the uncertainty of the background calculation to be 30%. 

TABLE II. Data, expected background, and expected stop 
signals after final cuts. Stop A scenario represents a small 
Amf_5 with mj {ma) = 100 (75) GeV/c^. Stop B scenario 
represents a large Arn^.^ with mi (mc) = 120 (60) GeV/c^. 



Selection 


Data 


Background 


Stop A 


Stop B 


Preselection 


176 


155.3 ± 50.2 


23.6 ± 8.9 


34.5 ± 13.0 


OS & E!t 


26 


28.7 ± 8.6 


12.9 ± 4.9 


25.1 ± 9.5 
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8.1 ± 2.4 


6.7 ± 2.5 


14.8 ± 5.6 


small Amf_ii 





1.5 ± 0.5 


5.7 ± 2.1 




large Amf_^ 





2.1 ± 0.5 




8.2 ± 3.1 




95% C.L. cross 


section limit: 


9.0 pb 


2.2 pb 



After establishing the selection cuts by using a "blind" 
analysis technique, we then apply the cuts to the pre- 
selection data. We observe zero events for both the 
small and the large Amf_^ sets of cuts, consistent with 
our background expectation of 1.52±0.47 and 2.07±0.46 
events. We used the frequentist method [23] with zero 
observed events, no background subtraction and a total 
uncertainty of 38% on the predicted signal to calculate 
a 95% confidence level upper limit of 4.01 stop events. 
Consequently, we exclude all stop-sneutrino mass combi- 
nations that would yield more than 4.01 events. Figure 2 
shows our result compared to LEP2 [6] and D0 [22]. The 
improvement in the small Amf_^ region is due to our in- 
creased signal to background ratio of 4:1. 

In conclusion, we have searched for stop pair produc- 
tion in 107 pb~^of data from pp collisions at ^/s = 1.8 
TeV collected by CDF. The observed dilepton, jet, and 
missing Et events are consistent with expectations from 
SM sources. Failing to find a signal of supersymme- 
try, we establish mass limits at 95% C.L.: we exclude 
stop masses up to = 135 GeV/c^ (at mp of 72- 
79 GeV/c^) and sneutrino masses up to 88.4 GeV/c^ (at 
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FIG. 2. Stop and sneutrino mass plane showing the CDF 
95% C.L. excluded region as hatched area. For the 3-body 
stop decay, t — > Wb, a 33.3% branching ratio to each of the 
three lepton flavors is used. 



mi of 126 GeV/c^). 
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